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Introduction 


The massless Schwinger model has been studied extensively because, although being exactly 
soluble, it has quite a rich structure that resembles features of more realistic theories. Among 
these are: a massive physical state is formed via the chiral anomaly 

[6j]) (this state is noninteracting in the massless Schwinger model). There are ”instanton-like” 
gauge field configurations present, and therefore the vacuum structure is nontrivial: the true 
vacuum is a superposition of all instanton sectors, and each possible superposition is labelled 
by a vacuum angle 9 (@. i, S, i, 0, i, i)- However the physics does not depend on 
the value of 9 in the massless case. 

The massive Schwinger model is no longer exactly soluble (Jl^], [11], [[HJ, 0]). But all of the 
nontrivial features of the massless model persist to hold, at least for small fermion mass. The 
massive state is now interacting, and its mass acquires corrections due to the fermion mass. 
Instantons and nontrivial vacuum are present, too, and, in addition, physical quantities now 
depend on the value of the vacuum angle 9. 

Here we will show how some physical quantities may be computed via a perturbative ex¬ 
pansion in the fermion mass, and we will arrive at some numerical results for the vacuum 
functional and fermion condensate. For this we use the wellknown exact path integral solution 
of the massless Schwinger model as a starting point. 


Perturbative computation 

The vacuum functional for the massive Schwinger model may be written as a sum of all instanton 
sectors (k... instanton number) 


oo 

Z(m,9) = J2 e iW Z k (m ) 

k=—o o 


(i) 


where 


Z k (m) = N 


DVDVDA^J** 




= Nj D*D*D0 k ^ ^ 11 / dx^{xi)^{x 


n =0 


i= 1' 


• exp 


dx 


+ -f3 k n 2 (5 k } 


( 2 ) 


(Ap = corresponding to Lorentz gauge) where we expanded the mass term. Therefore the 

vacuum functional is related to the computation of fermionic VEV of the massless Schwinger 
model (and in addition some space time integrations thereof). These fermionic VEV of the 
massless Schwinger model may be computed from the generating functional 


Zk[v,v] 



k—1 

Df3lUvK)«V> 

*o=0 


. e i / dxd vn(x)GP( x ,y)v(y) e \ f d.x/3 D/3 


( 3 ) 
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where the fermions have already been integrated out. are the fermionic zero modes and G 13 
is the exact fermionic Green’s function in the external gauge field /?, and D is the operator of the 
effective gauge field action, with Green’s function G (D^(x)(Dq(x)) is the massive (massless) 
scalar propagator) 

G(x) = 7 t(D im (x) - D q {x)) 

D n(x) = --^ K o(n\x\) , D 0 (x) = ^lnx 2 (4) 

(Kq ... McDonald function). From this all fermionic VEV may be computed (for details see e.g. 

®, i). 

It is useful to rewrite the fermionic scalar bilinears in terms of chiral bilinears like 


S(x) := 'F(x)4'(x) = 'F(x)P + 4 , (x) + 'F(x)P_'F(x) =: S+(x) + S-(x), 


(5) 


because for a VEV of products of S± only a definite instanton sector contributes. E.g. the 
fermion condensate is 


<S(*)> = (s+M)* -1 + (S-( x)) 


1 = h 2G «» = —a =: E 

7T 27T 


( 6 ) 


where 7 is the Euler constant and /i the Schwinger mass /j? = ^ 
For higher VEV of and n_5_ one finds 


(S Hl (x 1 ) • • • S Hn (x n )) = exp[^(-) CTi<T M 7 r - xj) 

i<j 


(7) 


where a t = ±1 for H, = ±. The sole contribution to this VEV stems from the instanton sector 
k = E, cr* = n+ - n_ (see e.g. [|, 0 ). 

Now we just have to insert this result into the perturbation expansion (2). For a first, 
rough approximation we may use the fact that the massive scalar propagator D^{x) vanishes 
exponentially for large argument. Therefore, when integrating over space time and expanding 
the exponential, all contributions from D l /t will be ignored for the moment, supposing that 
the space time volume V is sufficiently large. In this case the integrations in (2) just produce 
factors of V. Further, when inserting (7) in (2) we have to sum over all possible distributions 
of n+ = n — ri- pluses and n_ minuses on n scalar densities S. This results in a factor ( " ). 
Therefore we find for the n-th order term 


(II / dxiS(xi)) 


Tl 

(§)V V 




i =1 ri— =0 

and for the normalized vacuum functional 

Z(m,d) 


n 

. Tl- , 


e i(n-2n-)6 = cos Qy 


7(0,0) 


exp (mYiV cos 9) 


( 8 ) 


(9) 


(we ignored terms like m n V n ^ 1 in this approximation compared to m n ~ 1 V n ~ 1 , therefore (9) 
is the first order result in to) . This result is wellknown, and its consequences for the vacuum 
structure and spectrum of the Dirac operator are discussed in great detail in [18|. 
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To obtain higher order results, we rewrite the exponential in (7) like e.g. 
ex P(5Z(—) a%aj4:71 D,,.(x t - Xj)J = exp^+4vrT» M (xi - x 2 )'j ■ exp(-47rD /t (x'i - x 3 )) • ... 

= (l + E(xi - x 2 j S ) ■ (l + F(x i - x 3 )) • ... (10) 

where 

E[x) = e 4nD ^ x) - 1 , F(x) = e ~ 4nD ^ x) - 1 

E = J d 2 xE(x) , F = J d 2 xF(x). (11) 

Both functions E(x), F(x) vanish exponientially for large argument. 

Using this notation and inserting (7) into the perturbation expansion (2) we obtain, order 
by order: 


n=l: 


n=2: 


n=3: 


n 


=4: 


^!/cte(e«+e 


-»> = “| y ( e « + e -«; 


m f ^ 2 f d Xl dx9 e 2ie e 47rD v( x i~ x 2) _|_ 2e~ 4nD ^( Xl ~ X2 ' > + e ~ 2id e 4nD ^( x i~ x 2) 

2 ! ^2 * J 


... = ) 3 [U 3 (e 3ie + 3e ie + 3e~ id + e" 3ie )+ 

V 2 (3Ee 3ie + 3 (E + 2 F)e ie + 3 (E + 2 F)e~ id + 3 Ee~ m ^ + 
U((3 E 2 + E x E x U)e 3i0 + 3(2UF + F 2 + E x F x F)e ie + 
3(2UF + F 2 + E x F x F)e~ id + (3U 2 + £x£x U)e _; ’ 


■ 6 + 4e -2 * 0 + e 


~ m )+ 

~ 2ie + 6 Ee~ 4ie ) 


+ 


( 12 ) 


= (f ) 2 [ y2 ( e2ie + 2 + e ~ 2W ) + ^(^e 2ie + 2F + Ee~ 2ie )\ (13) 


(14) 


... = ^(f) 4 [r 4 (e“+ 4 e« + 

U 3 (6£e 4i0 + 12(U + F)e 2 * 0 + 12(U + 2F) + 12 (E + F)e‘ 

U 2 ((15U 2 +4BxEx E)e iie + 4(3 E 2 + 9 EF + 3F 2 + Ex£x£ + 3£xFx F)e 2ie + 
+6{E 2 + 8 EF + F 2 + AE x F x F) + 4(3 E 2 + 9 EF + 3F 2 + ExExE + 3ExFx F)e~ 2ie + 

+ (15£ 2 + 4£xfix£)e“ 4tf )+...] (15) 


where e.g. 


E x E x E = J dyidy 2 E(yi)E(y 1 + y 2 )E(y 2 ) 


(16) 
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and we displayed the result up to the accuracy we need. Observe that the result is not ob¬ 
tained by just expanding polynomials like (1 + E(xi) 1 , because e.g. a third power in E(x{) may 
contribute to V n ~ 3 E 3 or to V n ~ 2 E x E x E. 

In a next step we rearrange the terms (12) - (15) in rising powers of V: 

[m|(e“ + e-«) + ^(§ f(Ee m +2 F + Ee~™) + 

(-) ((3F 2 + ExEx E)(e 3ie + e ~ 3ie ) + 3(2 EF + F 2 + E x F x F){e ie + e ~ id )) +...] + 

[m 2 (|) 2 (e 2 « + 2 + e - 2 «) + m 3 (|) 3 (E(e 3 " + e - 3 ") + (E + 2E)(e" + e- 3 )) + 

(f ) 4 (< 15e2 + 4E * E « E)(e iie + e-")+ 

4(3F 2 + 9FF + 3F 2 + Ex£xB + 3f;xFx F)(e 2id + e~ 2id )+ 

6(F 2 + 8 EF + F 2 + 4F x F x F)) + ... + 

[ m3 (f ) 3 ( e3 * 0 + 3 ^ + 3e_ie + e“ 3i0 )+ 

^ (!)' 4 (3 E(e iie + e~ 4id ) + 6 (E + F) (e 2ie + e~ 2ie ) + (F + 2 F)) +...] + 

^ [m 4 (|) 4 (e 4ie + 4e 2 * 0 + 6 + 4e~ 2W + e~ 4ie ) + 

V V 2 o V 3 n 

= ■ l\ a + ^r a + "3!” a + ( 17 ) 

We find that the coefficient of -^r is the n-th power of some specific function a where a does not 
depend on V any more. This feature could have been expected, because the fermion condensate 
should not depend on the volume V. With this function a, 

^) 777 /^ / ^ ] \ 2 

a(m, 9) = m—2 cos & + J (2 E cos 29 + 2F)+ 

(-)' ((3F 2 + F x F x F) 2 cos 39 + 3(2 EF + F 2 + E x E x E)2 cos + ... (18) 

the normalized vacuum functional may be written like 

Z(m, 9) _ ^va{m,e) 

Z( 0,0) 

From this it is very easy to compute the fermion condensate 

(’F’FHm, 9) = — ——In Z(m,9) = —— aim, 9) 

V am, dm 

(M)(m,0) = S cos 9 + ^-£ 2 (F cos 29 + F)+ 


( 20 ) 
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^-£ 3 ((3 E 2 + E x E x E) cos 39 + 3(2 EF + F 2 + E x F x F) cos + ... (21) 

For a numerical evaluation of order m we have to compute the coefficients E and F. First, both 
E and F are proportional to For /t = 1, E is 

E = J d 2 xE(x) = J d 2 x(e~ 2Ko{lxl) - 1) 

poo 

= 27r / drr(e~ 2Ko(r) - 1) = -8.9139 (22) 

Jo 

E(x) is well behaving (-E(O) = —1), so the numerical integration is straight forward. F(x ) is 
singular at x = 0, F(x) ~ -4 for x —> 0, but this singularity can easily be understood and 
removed in a unique way. Indeed, this singularity is just the free fermion singularity, as can be 
seen by rewriting F{x) like 


Yf(F(x) + 1) = G?(*)e te < G < l X- G <-'» '*-!=? G§(z) = 
This singularity may be isolated by a partial integration: 

roo fa 


F = J d 2 x(e 2K °^ x ^ — i) = 2 7r ^ ( e 2ft -°( r ) + 2 

= 27r[lnr(e 2Ko(r)+21nr -r 2 )l° 0 

J €—>0 

poo I 

+ 2ir / dr2 In r((IF i (r) - _) e 2 ^o(r)+2 lnr + r) 
JO r 


lnr _ r “2\ 


(23) 


(24) 


(Kq = —K\). Observe that the first term precisely leads to the free field singularity at the lower 
boundary (and vanishes at the upper boundary). So the second term is the unique and finite 
result we are looking for. The numerical integration gives 


F = 9.7384 


(25) 


With these results we find for the fermion condensate 


(ff)(m,0) 


£cos0 H-^-(—8.9139 cos 26 + 9.7384) + o(m 2 ) 

2 fi z 


(26) 


or, using (6), 


-(W)(m,d) 

[1 


p T 777 p 777 ^ 

— cos e + - — (E cos 20 + F) + o(—«-) 

Z7T /i 87T Z /i^ 

9 

in Tfi 

0.2835 cos 6 + —(0.7825 - 0.7163 cos 26>) + o(^-). 

[i [i z 


( 27 ) 


This is our final result. Observe that the correction is minimal for 6 = 0. 



6 


Summary 

We have derived a systematic procedure of computing the vacuum functional and fermion 
condensate of the massive Schwinger model in a mass perturbation theory, order by order. For 
the first nontrivial order we even displayed numerical results. For higher orders the numerical 
calculations are more involved (because of the occurrence of linked functions in the integrations, 
like E x E x E), but in principle they could be performed, e.g. in order to compare the result 
to the results of other approaches. 

The numerical computations in this article were done with Mathematica 2.2. 
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